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FOREWORD 


The objective of the four-dimensional atmospheric modeling task 
is to provide world-wide profiles of pressure, temperature, density, and 
moisture from the surface to 25 km altitude. The model gives mean 
monthly values and daily variations of the four parameters for any input 
of latitude, longitude, and month. This model can then serve as input for 
attenuation models that predict the degree of atmospheric attenuation 
likely to be encountered by satellite or air-borne electromagnetic sensors 
engaged in earth resources observations. It can also be used as mean 
model atmospheres in trajectory and vehicle heating analyses. 

This report gives a comparison study between the 4-D model versus 
standard model atmospheres. Also radiometric computations using these 
two models are given and discussed. 

The 4-D computer program can be obtained, upon request, from 
the Aerospace Environment Division, NASA, Marshall Space Flight Center, 
Alabama 35812. 

Associated work has been and is still being done in this field. For 
example, work reported on, thus far, includes the world-wide cloud cover 
model (NASA CR-61226, NASA CR-61345, and NASA CR-61389, and the 
interaction model involving microwave energy and atmospheric variables 
(NASA CR -61348). The four-dimensional atmospheric model is being 
improved to ultimately present pressure, temperature, density, moisture, 
and cloud cover as one attenuation model for earth resources problem. 
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1. INTRODUCTION 


The Aerospace Environment Division (Aero-Astrodynamics Laboratory) of the 
Marshall Space Flight Center has sponsored the development of a global atmos- 

* pheric model, known as the MSFC 4-D Atmospheric Model. This Model consists 
of 45 "homogeneous" moisture regions; for each region, pressure, temperature, 

. density, and water vapor are given at kilometer levels up to 25 km. Previous 

studies (Spiegler and Greaves, 1971; Spiegler and Fowler, 1972) have processed 
the necessary global data and refined the atmospheric statistics; this study 
has attempted to analyze the meteorological and radiometric differences in- 
dicated by the region values. 

Regions in the 4-D atmospheric model were defined on the basis of atmo- 
spheric moisture; each region is homogeneous with respect to its moisture 
mean and seasonal variability. As water vapor amounts are dependent on 
atmospheric temperature and reflect general climatic patterns, these moisture 
regions are also presumed relatively homogeneous with respect to the tem- 
perature, pressure and density. The atmospheric statistics given for these 
regions consist of monthly means and variances for each meteorological parameter; 
these values are derived from similar statistics given for 3470 points span- 
ning the globe. The grid point values form a more detailed version of the 4-D 
atmospheric model; a brief description of these models, and of the data used 
to derive them, is given in Table 1-1. The U.S. Standard Atmosphere Supplement , 

1966, defines representative atmospheric profiles on the basis of "latitudinal 
homogeneity"; the first part of this study was devoted to a comparison of 
these two atmospheric models, and to an analysis of the meteorological aspects 
of this "latitudinal homogeneity." 

The second part of this study investigated the differences in the two 
atmospheric models in their simulation of remote sensing measurements. The 
increased use of remote sensing data by individuals not skilled in the in- 

* terpretation of atmospheric effects necessitates tables which provide easily 
accessible correction values. To date, such tables have been based on the 
Standard Atmospheres; in this study, the value of tables computed from 
more detailed atmospheric data was investigated. This analysis paralleled 
the study of atmospheric latitudinal homogeneity, and was summarized by two 
simulated orbital paths. 
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To perform the above studies, certain minor investigations were neces- 
sary. First, conversion was needed from the 4-D atmospheric model's 
moisture parameter of absolute humidity, to the more common parameters of 
mixing ratio and dew point. Appendix A summarizes the technique used. 
Second, no new radiometric model was developed for this study; rather, the 
programs developed by Dr. David Anding were adapted to analyze the Standard 
Atmosphere and the 4-D atmospheric model. Appendix B reviews these mo- 
difications. The third appendix, Appendix C, presents an evaluation of 
the effects of clouds on the moisture statistics to permit an estimation 
of the applicability of these models to clear day situations. 

TABLE 1-1 

SUMMARY OF THE 4-D ATMOSPHERIC MODELS 

INPUT DATA TO 4-D MODELS 


Parameters 

Form 

Levels 

Grid 

Number of Points 

Height 
Temperature 
Dew Point 

Daily 

Values 

(Eight Years) 

1000 - 300mb 
1000 - 3Q0mb 
8S0 - SOOmb 

The National 
Meteorological 
Center (NMC) Grid 

1977 

Height 

Temperature 

Daily 

Values 

CSeven Years) 

100 - lOmb 
100 - lOmb 

NMC 

Grid 

1977 

Height 
Temperature 
Dew Point 

Monthly 
Means 5 
Variances 

Surface - lOOmb 
Surface - lOOmb 
Surface - SOOmb 

5° Latitude - 
Longitude Grid 
Northern Hemisphere 

1368 

Height 
Temperature 
Dew Point 

Monthly 
Means S 
Variances 

Surface - lOOmb 
Surface - lOOmb 
Surface - SOOmb 

5° Latitude - 
Longitude Grid 
Southern Hemisphere 

1368 

Height 
Temperature 
Dew Point 

Daily 

Variances 

850 - lOOmb 
850 - lOOmb 
850 - SOOmb 

5° Latitude - 
Longitude Grid 
0° - 15° North 

288 


FINAL 4-D MODELS 


MODEL A: Global Grid Point Data 

Parameters 

Form 

Levels 

Grid 

Number of Points 

Pressure 
Temperature 
Absolute Humidity 
Density 

Means 6 
Dai ly 

Variances for 
Each Month 

Surface - 25 km 

NMC 

Grid 

1977 

Pressure 
Temperature 
Absolute Humidity 
Density 

Means § 
Daily 

Variances for 
Each Month 

Surface - 25 km 

5° Latitude - 
Longitude 
5° - 90° South 

288 

* MODEL B: Global Moisture Regions 


Pressure 
Temperature 
Absolute Humidity 
Density 

Region Means 
5 Daily 
Variances for 
Each Month 

Surface - 25 km 

"Homogeneous” 
Moisture Regions 
Entire Globe 

45 
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2. REGIONAL DIFFERENCE IN THE 4-D MODELS 

Previous studies have established the homogeneity within each region 
of the 4-D atmospheric model. This section will investigate the homogeneity 
^ within each latitudinal belt of the U. S. Standard Atmosphere Model. The 

analysis will be based primarily upon regional differences in precipitable 
water amounts, but will also investigate differences in the vertical structure 
of temperature and dew point. 

2.1 Region Definition 

It is customary to describe the atmosphere over a given location by its 
latitude. Thus Washington, D.C. has a mid-latitude atmosphere, Trinidad a 
tropical one, and Thule, Greenland is considered subarctic. These descrip- 
tions assume that each latitudinal belt is homogeneous with respect to 
temperature, moisture and pressure, and that a knowledge of the atmospheric 
parameters for any point on a given latitude implies a knowledge of condi- 
tions for all points on that latitude. As this homogeneity is approximately 
true, the Handbook of Geophysics provides atmospheric profiles representa- 
tive of given latitudes and seasons. These will be referred to below as 
the standard atmospheres. 

In reality, latitudinal homogeneity is affected by the underlying 
surface. This surface is not homogeneous since oceans and continents are 
irregularly spaced and the topography of the continents varies markedly. 
Furthermore, the vast currents of the atmosphere and the oceans are not 
entirely zonal in character. To permit a study of the effects of these 
factors on atmospheric moisture, homogeneous water vapor regions (Table 
2-1 and Figure 2-1) were defined based on the total amount of precipitable 
j water (Figures 2-2 and 2-3) and its seasonal variability (Spiegler and 

Greaves, 1971). 

As the table and figures show, there are 45 regions in all, in contrast 
to the five given by the U.S. Standard Atmosphere Supplement , 1966. 

Some of these regions are zonal, such as the maritime ones, but others are 
meridional (reflecting major mountain ranges), and others are applicable only 
to small areas of the globe such as the Sahara or the Australian Desert. 

Each region differs from all other regions, either in its mean value of 
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TABLE 2-1 

HOMOGENEOUS MOISTURE REGIONS 


Region 

Number 

bbi 

IBB) 

r i .tli ; 1 ily A< r> ss Rigi'.i. 

I . if < 1 1< in 

Remark s 

1 

Flxtrcmoly low 

Small 

Small 

Arrl - 


2 

Very Urn 

Small to inoler.itc 

sm. ill, except moderate in 

S'di.i |., t it ..i',l Arctic - 
i "id me lit 


» 

I ow 

Si * i « t 1 1 

Small, except moderate in 

Snh.irrl ii 

V «• r v low lorn Mu r*- late fall to 
early spring 

4 

I.I1W 

Murler ite to largo 

Sm .ill, except large in 
summer 

1 'ol.tr - < , . nt i m<nl 

Moderate moisture summer 

S 


Mod. !,. r l; ;• 

Small, except moderate m 

N'il>p.,|.,r 

Very low musture winter 

i. 


Sir. HI 

S'" 'll. '-v t opt I. .rye in 

1 'ol.i r <■. 


7 

I.nw to moderate 

Mod, rale 

Small I. ill, winter and 
early spring, largo in 

High Ol id - 1 .if ilnrli' - iindl- 
nioit 

Summer maximum higher, and 
Regions H, “ and 10 

8 

Low to moderate 

Small to moderate 

Sm.itt In moderate 

t iigt, II 1 1 1 - 1 ,f ifuilr IH IMII 


9 

Lnw fu moderate 

Small to moderate 

small, cm cpI moderate 
slimmer 

High mi.| - l.itii ode mean 

Winter maximum lower than 
Region K 

10 

Low to moderate 

Small 



High mid - 1 .itii ndt- - ocean 

Moisture higher than Regions H >. ') 

1 1 

Moderate 

Small to moderate 

Moderate! oxrcpt la v\*v m 
jiiimmi'r 

Mi d-l.ititndi - in can 


12 

Moderate 

Small to moderate 

Sm al 1 to moderate 

low nwH - l.ih 

nebcin* 

Higher winter moisture tlian 
Region 11 

n 

Moderate 

Large 

Moderate, except small 
n winter and large early 
■ all 

Mid- 1 .tlitudc continent 


14 

Moderate 

Large 

Moderate, except large In 
(all 

Lower mid-latitude conti- 
nent 

Somewhat higher average than 
Regions 1 1- M. 

15 

Moderate 

Moderate to large 

Moderate, except largo in 

fall 

Lower mid- latitude 
oceans 

V cry high summer and early fall 
moisture; winter moisture higner 
than Region 14 (lowest 2 km) 

16 

Moderate 

Small to moderate 

Moderate 

Snhtropic.il eastern por- 
tion ol oceans 


17 

Moderate to high 

Moderate 

Moderate to large 

Subtropics including 
western 2»'t of m oans 


18 

Moderate to high 

Small 

Small, oxcept email to 
moderate in summer and 
fall 

Eastern subtropical and 
tropical oceans 

W inter moisture higher than 
Region 17 

10 

High 

Small to moderate 

Moderate, exrcpt large 
in aprmg and fall 

Subtropical - tro pi cal 
horde r 


20 

High 

Small 

Small, except moderate 
in summer and tall 

Tropical eastern portion 
of oceans 


21 

Very high 

Small to moderate 

Moderate, except large in 
■ pring 

Tropical 


22 

I 

Very high 

Small 

Moderate 

Tropical oceans except 
near Ec) . in E. oceans 


23 

Moderate 

Moderate 

Moderate, except small 
m spring 

Tropical interior conti- 
nent region 


24 

Extremely high 

Very small 

Small 

Equatorial oceans tor 
parts of globe 


25 

Moderate to high 

Small to moderate 

Small, except moderate 
in winter and spring 

Equatorial continent 

Lower winter moisture than Region 
24 

26 

Low 

Small to moderate 

Small to moderate 

M.uor mountain ranges 
average ulev. •*> .!-3 km 


27 

Very low 

Small 

Small to moderate 

Major mountain ranges 
average clev. -*.4-7 km 


28 

Low 

Small 

Small 

border of mountains mid 
and high latitudes 


29 

Low to moderate 

Moderate to large 

Moderate, except large in 
aumin er 

border ol mountains sub- 
tropics and tropics 


30 

Low 

Sm all 

Small, except moderate 
in iummer 

Deserts 


11 

Low 

Small to moderate 

Small 

border ol desert 


32 

Low to moderate 

Moderate 

Small to moderate 

Border oi desert 


33 

i-ow to moderate 

Moderate 

Moderate, except large in 
fall 

border of desert 
(tropica) and equatorial! 

Average moisture similar to Region 

14 

Moderate 

Very large 

Very large, except email 
in summer 

Mid - 1 at i tude - E, coast 

Asia 


15 

Moderate to high 

Small to moderate 

Moderate to large, except 
small in fall 

Tropii « - near continent- 

Similar to Region 19, except lower 
moisture in spring 

36 

Moderate to high 

Large 

Moderate, except small in 
summer 

India and Indian Ocean 

Extremely high late spring and 
summer moisture 
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TABLE 2-1 Continued 


Region 

Avernge Annual 
Moisture 

Seasonal 

Change 

Variability 
Across Region 

Location 

Remarks 

37(1) 

Very low 

. 

Very small 

Small 

Antarctic continent 

Higher annual average mois- 
ture than arctic N. H. 

(Region l) except lower fat 
summer. 

38(2) 

Very low 

Very small 

Small 

Northern Antarctic 
Ocean 

Higher annual average mois- 
ture than subarctic and 
arctic continent N. H. 
(Region 2) winter and spring; 
lower in summer. 

39(6) 

Low 

Very small 

Small 

S. H. polar ocean 

Higher annual average mois- 
ture than N. H. polar oceans 
(Region 6) except in summer. 

40(7) 

Low -to -mode rate 

Very small 

Very small 

Southern one-third 
South America 

Lower annual average mois- 
ture than high midlatitude 
continent (Region 7), except 
in winter. 

41(15) 

Moderate 

Small -to-mode rate 

Small 

S. H. lower midlati- 
tude oceans. 


42(26) 

Low 1 - to -moderate 

Small -to -moderate 

Small -to -moderate 

S. H. mountain ranges. 


43(30) 

mm 

Very small 

Small 

Australian Desert 

Moisture lower than average 
of N. H. deserts (Region 30) 
in summer. 

44(31) 

Low 

Small 

Small 

Border of Australian 
Desert. 

Temperature much lower than 
N. H. border of desert (Re- 
gion 3 1). 

45(33) 

Low -to -moderate 

Moderate 

! 

Moderate except 
large in fall. j 

Tropical border of 
Australian Desert. 

Lower moisture than N. H. 
tropical and equatorial desert 
borders (Region 33). 


































































moisture or in its variability. While such detailed differences are re- 
flections of actual climatic differences, it is not immediately clear that 
such precise region definitions are required for the intended application. 

There is, therefore, a need to investigate the degree of variability in 
the moisture data for all regions falling within the latitudinal belt of 
one of the standard atmospheres. Such a study was performed using the 
January and July data from the 4-D models; the following paragraphs present 
the results. 

2.2 Comparative Analysis of Moisture Data 

For this analysis, the nine atmospheres applicable to 15°N, 30°N, 45°N, 60°N, 
and 75 °N were selected for January and July* from the Handbook of Geophysics. 

The 15°N atmosphere is assumed applicable to the latitudes from 0°-22.5°, 
the 30°N, to 22.5°-37.5°, the 45°N, to 37.5°-50.25°, the 60°N, to 50.25°-60.75°, 
and the 75°N from 60.75° to the pole. These latitudinal belts are considered 
appropriate for both hemispheres although the seasonal reversal between 
hemispheres necessitates the use of January data for the Southern Hemisphere . in 
July, and the use of July data for January. These atmospheres are shown in 
Figures 2-4 to 2-13 as temperature-dew point profiles; with them are presented 
similar profiles of the region parameters to allow a comparison of the data at 
various latitudes. Total precipitable water was also computed for the standard 
atmospheres, and for all regions, and the resulting values were plotted for each 
latitudinal belt as shown in Figures 2-14 to 2-23. 

2.3 Discussion 

2.3.1 Latitudinal Belt from 60.75° to the Pole 

I Examination of the standard atmospheres and regions applicable to 

latitudes greater than 60.75° (Figures 2-4, 2-5, 2-14, 2-15) makes several 

- facts immediately clear. The most obvious is the marked difference between 

the two hemispheres. The Antarctic winter is so much more moist than the 
Arctic winter, and the Antarctic summer so much drier, that the standard 
atmospheres for 75°N (or any other atmospheres applicable to the Northern 
Hemisphere) are not representative for either Region 37 or 38. 

*0nly an annual profile is given for 15°N because of the limited seasonal 
variability. 
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Figure 2-2 Global Precipitable Water (g cm ) January 








Figure 2-3 Global Precipitable Water (g cm ) July 






Figure 2-4 Temperature and Dew Point Profiles 60.75° - 90 N January, 
60.75°- 90°S July 






Pressure (mb) 



- 40 ° - 30 ° - 20 ° - 10 ° 0 ° 10 ° 

Temperature (°C) 


Figure 2-5 Temperature and Dew Point Profiles 60.75° -90°N July, 
60.75° -90°S January 
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W JULY 



Temperature (°C) 

Figure 2-7 Temperature and Dew Point Profiles 50.25° -60.75°N July, 
50.25° -60.75°S January 
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Figure 2-10(a) Temperature and Dew Point Profiles 22.5 -37.5 N January 
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Figure 2- 11 (a) Temperature and Dew Point Profiles 22.5 -37.5 N July 

22.5° -37.5°S January 
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Figure 2-llCb) 
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Figure 2-12 Temperature and Dew Point Profiles 0 -22.5 N January, 

0° -22 . 5°S July 
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Figure 2-13 Temperature and Dew Point Profiles 0° -22.5°N July 
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Figure 2-14 Precipitable Water 60.75 
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Figure 2-16 Precipitable Water 50.25 -60.75 N 
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Figure 2-19 Precipitable Water 37.5 -50.25 S 
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Figure 2-20 Precipitable Water 22.5° -37.5°N 
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Figure 2-23 Precipitable Water 0° -22.5°S 


3 


A temperature inversion appears in the January standard atmosphere 
corresponding well to data recorded in the high latitudes. The inversion 
is caused by the radiational cooling of the snow-covered surface below, 
but averaging across time and space results in its disappearance in the 
region data. The dew point inversion, however, is retained in Regions 1 

and 2 (the Arctic regions), and the standard atmosphere's precipitable 

-2 -2 
water (0.19 g cm ) falls about halfway between Region l's 0.15 g cm and 

Region 2*s 0.25 g cm “. The July standard atmosphere's profile is unlike 

the region profiles for both temperature and dew point, yet the temperatures, 

and precipitable water, again fall about halfway between those of the two 

regions . 

Region 3, the subarctic ocean, is also located predominantly north 
of 60.75°N. It is, however, much warmer and more moist (Figures 2-6, 2-7 
and 2-8 than either of the Arctic regions or the standard atmosphere at 

75°N. Thus the latitudes between 60.75° and the poles would be best re- 
presented by the region data, although use of the 75 °N standard atmosphere 
for the Arctic regions would not cause large errors. 

2.3.2 Latitudinal Belt Between 50.25° and 60.75° 

The region between 50.25°N and 60.75°N is a transition zone from the 
high mid-latitudes to the subarctic ocean; the most appropriate standard 
atmospheres would be those at 60°N. The Southern Hemisphere at these 
latitudes is almost entirely ocean, implying that the variation around the 
globe in this zone is limited. However, the same hemispheric moisture 
differences mentioned for the Antarctic are evident in Region 39 (Figures 
2-6, 2-7, 2-16, 2-17), making the standard atmosphere for 60°N clearly in- 
appropriate. 

The Northern Hemisphere presents a more diverse picture with preci- 

_2 

pitable water values in this belt ranging from 0.25 to 1.1 g cm in the 

_2 

winter and from 1.75 to 3.22 g cm in the summer. The January standard 
atmosphere is most representative of land conditions at this latitude, with 

the temperature and dew point values corresponding well to Region 4 (polar- 

• -2 

continent) and the precipitable water value of 0.42 g cm matching Region 

_2 

4' s 0.41 g cm . Although slightly lower in temperature and water vapor 
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content, it is also quite similar to Region 5. However, the standard atmo- 
sphere greatly underestimates the temperatures and moistures found in sub- 
arctic and polar oceans (Regions 3 and 6) . The standard atmosphere pro- 
file also shows a slight inversion at 1 km; this is a characteristic of 
land or ice-capped regions only. 

The July standard atmosphere resembles land conditions only in the 
shape of its temperature profile since its temperatures run 5°C colder 
than those above Regions 4 and 5 and its dew point values are almost as 
low below 700 mb. In this season, the temperature values are much more 
similar to those in the subarctic regions, and are not unlike those of 
Region 6 except at the surface. The dew point profile, however, has a 
much smaller lapse rate than is found for any related region; it also 
underestimates the moisture for all but the subarctic regions (2 and 3) . 
Clearly, the variability around this latitudinal belt recommends use of 
regional data, rather than one standard atmosphere. 

2.3.3 Latitudinal Belt Between 37.5° and 50.25° 

The standard atmosphere for 45° is again most closely allied to land 
(Figures 2-8, 2-9, 2-18, 2-19). Like the atmospheres for higher latitudes, 
it underestimates the winter moisture, and overestimates the summer moisture 
throughout the Southern Hemisphere regions from 37.5° to 50.25°S, and should 
never be used there. However, its temperatures and summer dew points well 
represent Region 7 (high mid- latitude continent) and Region 28 (border of 
mid- and high-latitude mountains), and the precipitable water, though a 
little low, would not provide a bad estimate of the values in those two 
regions. The same cannot be said of the high mid-latitude oceanic regions 

8 and 9. In these areas the winter temperatures run higher than the standard 

atmosphere's and the summer temperatures run lower. Likewise the precipi- 
table water values are much higher in January and much lower in July. The 
mid- latitude ocean region. Region 11, represents a cross between these land 
and ocean climates. In the winter, its climate resembles those of Regions 
8 and 9 with somewhat higher temper atures and moisture; in the summer, its 
temperatures coincide with the 45°N standard atmosphere’s, and its precipi- 
table water almost equals that of Region 7 (3.0 g cm vs 3.2 g cm ), This 
drastic change in climate is not surprising; in January this region 
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is under the influence of the Icelandic low (as are Regions 8 and 9) , in July 
it is subjected to warm, moist air drawn north by the circulation of the 
Bermuda high (as is Region 7) . But it underlines not only the variability 
of moisture within this latitudinal belt, but also the seasonal variability 
of that zonal variability. 

2.3.4 Latitudinal Belt Between 22.5° to 37.5° 

The latitudinal belts from 22.5° to 37.5°N and from 22.5° to 37.5°S are 
characterized by great variability in geography. In the Northern Hemisphere 
these latitudes include the Sierra Madres, the subtropical oceans, the 
Sahara Desert, and the Himalayas. In the Southern Hemisphere these 
latitudes include the southern half of Australia, the mountain ranges 
of Chile, the plains of Argentina and the Union of South Africa. Almost 
every type of climate is represented from the arid, to the tropical wet, 
to the highland, to the continental temperate. Clearly no one atmosphere, 
such as the standard atmosphere for 30°N, could be representative of the 
entire zone (Figures 2-10, 2-11, 2-20, 2-21). 

Precipitable water in the Northern Hemisphere ranges from 0.4 to 3.4 

-2 -2 
g cm in January, and from 2.15 to 6.6 g cm in July. The Southern 

-2 

Hemisphere is slightly less variable; its values range from 1.1 to 3.4 g cm 

-2 

in July, and from 2.0 to 5.35 g cm in January. The limited amount of land in 

this hemisphere is probably responsible for the somewhat smaller range of values. 

-2 

The winter value of 2.1 g cm , given by the standard atmosphere, is obviously 
a midpoint value; however, it is really only representative of Regions 10, 

12, 15 and 41, the lower mid-latitude oceans, and of Region 33, the border 
of Northern Hemisphere deserts. The standard atmosphere underestimates the 
moisture above the subtropical and tropical oceans, and above the continental 
tropics and subtropics (Regions 29, 35 and 36). On the other hand, it 
overestimates the moisture over all other land regions located in that 
latitudinal belt. 

_2 

The 30°N standard atmosphere for July contains 4.1 g cm of precipitable 
water, placing it with the most moist of the regions between 22.5° and 37.5°. 
Again it approximates the moisture in Regions 15 and 33, but it greatly over- 
estimates the moisture in Regions 10, 12, and 41. The value is also 
representative of the moisture in Region 14 (lower mid- latitude continent) 
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and in Region 26 (major mountain ranges), both regions that the January's 

standard atmosphere overestimates, and it matches well with the value 

for Region 17 (subtropic ocean) which the winter value underestimates. The 

only regions whose moisture values exceed the standard atmosphere's are 

Regions 19, 21, 29, 34, 35, and 36. Region 36 is composed of India and the 

Indian Ocean; in July it is in the midst of its monsoon season and has the 
• . -2 

highest precipitable water (6.6 g cm ) recorded for any season at any location. 

Nor surprisingly, the driest atmospheres in the subtropics are associated 
with the desert borders and the mountain ranges. 

The temperature profile given by the standard atmosphere for January 
seems to be a cross between that above the "eastern low mid- latitude ocean" 

(Region 12) and that over the deserts (Region 30) . Below 2 km it runs 
warmer than the temperature over the lower mid- latitude continent (Regions 14) 
yet it is consistently colder than the atmosphere over the subtropical oceans 
(Regions 16 and 17) and over the subtropical mountain borders (Region 29) . 

In the Slimmer the standard atmosphere's temperatures match those of Region 14 
near the surface, those of Regions 16 and 17 between 2 and 5 km, and are 
too warm above that level for any values but those of Region 29. Below 
2 km the standard atmosphere temper at ures exceed those over all the oceanic 
regions and fall short of those recorded in the deserts. 

2.3.5 Tropical Regions 

The tropical latitudes from 22.5°N to 22.5°S (Figures 2-12, 2-13, 2-22, 

2-23) (the Tropic of Cancer to the Tropic of Capricorn) share many moisture 
regions with the subtropics. In general, however, the latitudes do not show ex- 
treme seasonal variabilities; a fact which caused only an annual standard atmosphere 
to be computed for 15°N. Except for a slight inversion between 2 and 2.5 km, 
the temperature profile is quite similar to that given for 30°N July; thus 
the above comments about its applicability to the deserts and mountain regions 
still hold. No inversion appears in any region profile, presumably because 
of the temporal and spatial averaging, but the standard atmosphere temperatures 
below 4 km agree well with those of Region 22 (the tropical oceans) for both 
seasons. The other tropical regions (20, 21, and 23) are colder in January 
near the surface; all four tropical regions are warmer above 5 km than the 
standard atmosphere would indicate. 
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The dew point profile is totally nonrepresentative, especially as it 

shows a sudden, sharp drop in moisture between 2 and 2.5 km that is not 

even suggested by the region data. The total precipitable water given by 

_2 

the standard atmosphere, 3.9 g cm , overestimates all winter region values 
except those of Regions 22, 24 and 25. On the other hand, it underestimates 
all the summer values except for those of the subtropical oceans (Regions 12 
and 16), the mountain ranges (Region 27), and the deserts (Regions 30, 31, 32, 
42, and 45). 

Despite the inapplicability of the standard atmosphere as the represen- 
tative of this latitudinal belt, the spatial variability of moisture over most 
of this zone is not large. The amount of land found in the Northern Hemisphere 
south of 22.5°N is quite limited since Africa is the only continent that does 
not lie mainly to the north; probably two-thirds of this area is ocean. 

Oceans, of course, dominate the Southern Hemisphere, and though this belt con- 
tains much of South America and a good piece of Africa, less than one-fourth 

of these latitudes is occupied by land. Thus the winter precipitable water 

-2 -2 

range of 1.05 to 4.8 g cm can be reduced to 3.0 to 4.8 g cm merely by con- 
centrating on the regions which are predominate at those latitudes. Region 19 

_ 2 

through 22, and Region 24. Likewise the summer range of 2.15 to 6.6 g cm 

_2 

can be, without much error, regarded as 4.25 to 5.25 g cm . 

2.4 Summary of Comparison of Moisture Regions with 
U. S. Standard Atmospheres 

Although the significance of the regional atmospheric difference can only 
be evaluated in light of particular applications, the above analyses can be 
summarized as follows: 

• The greatest spatial variability within any latitudinal belt occurs 

between 22.5 # N and 37.5°N where 36 of the 45 regions are represented. This 

_2 

area has July values of precipitable water ranging from 2.15 to 6.6 g cm 
and many areas are seriously misrepresented by the 30°N standard atmosphere. 

• North of 37.5°N, the standard atmosphere is much more representative 
of land conditions than of ocean conditions. Generally, continental areas 
can be approximated by the standard atmospheres; the ocean regions are 
warmer and wetter in the winter, colder and drier in the summer. 
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• The U. S. Standard Atmospheres are not applicable to Southern 
Hemisphere latitudes south of 40°S. 

• The tropics show the least seasonal variability; however, the 
annual value given for 15 °N does not adequately represent these climates 
It overestimates the winter moisture in most of the regions and usually 
underestimates the summer moisture. Separate atmospheres for January 
and July would better represent this area. 
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3. THE 4-D MODELS AND RADIATION TRANSFER 


In this section, the use of atmospheric models in remote sensing 
studies is discussed in terms of a simple model of the transfer of infra- 
red radiation. This is followed by a short description of the computation 
technique actually used to convert the atmospheric data into the radiometric 
parameters that will be discussed in Section 4. 

3.1 Elements of Radiation Transfer 

In remote sensing by electromagnetic means, measurements are made 
of the radiation emitted, reflected, or transmitted by a target using a 
receiver located at some distance from the target. Properties of the 
target are then inferred from the radiometric measurements using known target 
signatures. Such signatures may include the intensity of the signal at one 
wavelength; it may include multispectral properties in which the relative 
intensities at a number of wavelengths are important; or it may extend into 
the angular domain in which the relative intensities measured at a number of 
angles provide the clue for target identification. The success of remote 
sensing techniques, whether to obtain an absolute measurement of a physical 
property of a target such as its temperature, or to simply identify the 
target, requires a consistent and known relationship between the measured 
signal and the properties of the target. 

With few exceptions, remote sensing measurements are made through the 
atmosphere. In the case of remote sensing of earth resources from an orbiting 
spacecraft, the sensors on the spacecraft measure the radiation from the target 
after the radiation has traveled through the entire atmosphere. This 
atmosphere, its gaseous and particulate constituents, are "radiometrically 
active". They absorb, emit and scatter radiation. In so doing, they 
interact with the radiation propagating from the target to the sensor. 

If this effect is so strong that inherent relationships between target and 
emitted (or reflected) radiance is totally destroyed or masked, then the ra- 
diation energy measured at the sensor cannot provide any information about the 
target. On the other hand, if the interaction is weak, the remotely 
measured radiation energy contains useful information concerning the target, 
and the measured signal can be used to infer information about the target. 

The degree of success in this inference is then dependent on the degree 
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to which the atmospheric effects on the measured signals is known and can 
be accounted for. 

More specifically the above discussion may be expressed in terms of a 
simplified equation applicable to the transfer of infrared radiation through 
an absorbing (and emitting) atmosphere as 

f ^.( 2 ) 

N x = E x (s) t x (A) N x [T(s)l + J N x [T(z)j — dz (3-1) 

0 


N x is the effective spectral radiance measured by a spacecraft borne sensor 
at wavelength X; E x (s) and T(s) are respectively the spectral emissivity and 
the temperature of the surface; and t x (A) is the spectral transmissivity of 
the entire atmospheric layer between the target surface and the sensor along 
the optical path. N x (T(s)J is then the blackbody spectral radiance at tem- 
perature [T(s)]. T(z) and t(z) are respectively the temperature and 
spectral transmissivity of the atmosphere at level z. 

Given the above definitions, Equation 3-1 may be simply interpreted 
as follows: the radiance measured at a satellite sensor consists of two 

parts. The first part is that which is emitted by the surface and 
attenuated by the atmosphere. This contribution is expressed by the first 
term in the right hand side of Equation (3-1) . The second contribution, 
expressed by the second terms in the equation, is made by the atmosphere. 

The total effects of the atmosphere on the measured radiance are through 
the terms of t x (A), t x (z) and T(z). It is in the evaluation of these 

terms that model atmospheres, such as the 4-D Models, are used. 

From Equation (3-1) it is clear that in radiation transfer studies 
the temperature profile of the model atmospheres is used implicitly through 
the blackbody Planck Radiance function. This same atmospheric parameter 
is also used in the evaluation of t (z) and therefore t (a) which is the 

A A 

integrated effect of the entire atmospheric path. Other parameters of 
the atmospheric models are also used in the evaluations of this term. 

Absorption or transmission of radiation by the atmosphere is spectral- 
ly dependent. Given the gaseous makeup of the atmosphere, at certain 
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wavelengths the atmosphere is approximately transparent while at other 
wavelengths, the atmosphere absorbs. The absorption by the different 
species of gases and the magnitude of the absorption is dependent on the 
amount of gas through which the radiation must propagate. In the infrared 
region, two gases are of primary importance. The first is water vapor 
and the second is carbon dioxide. It is clear, therefore, that to obtain 
values of t x (z), and of ^(A), it is necessary to use appropriate values 
of the mixing ratio profile W(z) or some other water vapor profile 
parameter. For a well-mixed gas such as carbon dioxide (or oxygen), only 
the pressure-height (or density-height) data from the appropriate atmospheric 
model are required to completely specify its mixing ratio profile (at least 
within the region of the atmosphere in which the gas is well-mixed). Yet 
a third atmospheric parameter is used in the evaluation of x^(z). This 
is the pressure profile, P(z). Due to pressure broadening effects on in- 
dividual absorption lines, the degree of atmospheric absorption is depen- 
dent on altitude not just because absorber concentrations generally de- 
crease with height, but also because pressure is a function of height. 

From the purely heuristic discussion given above, it should be clear 
that simulation of expected radiometric signals is possible, and often used, 
to evaluate potential problems in data interpretation of remote sensing 
measurements. More importantly, such simulations require representative 
values of detailed structures of the atmosphere such as might be obtained 
through standard atmosphere models or from the more detailed 4-D model 
atmospheres discussed in Section 2. 

3.2 Spectral Region Selection 

This study has addressed itself to the evaluation of the use of the 
4-D models in the infrared. The infrared region was selected because of 
its extensive use in quantitative measurements of surface and atmospheric 
properties in aircraft and satellite programs. More importantly, this 
region has significant water vapor absorption bands. Water vapor is the 
most variable of the atmospheric parameters. It is in fact the variabi- 
lity, spatial and temporal, of this parameter which has resulted in the 
definition of the large numbers of "homogeneous" regions in the 4-D atmospheres. 
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By confining the analysis to the infrared, it is hoped that validity or 
usefulness of the detailed 4-D models can be evaluated with a minimum number 
of computations . 

Computations were made in the following four spectral intervals: 

(1) 3-4.5 ym; (2) 6-8 ym; (3) 9-13 ym; and (4) 13-18 ym. 

The first and third of the spectral intervals are located in the so- 
called atmospheric windows where atmospheric effects are minimal, and have, 
therefore, been used for surface measurements, including surface temperatures. 
However, residual absorption by water vapor and carbon dioxide is still 
present. The first wavelength interval corresponds to the spectral region 
of the High Resolution Infrared Radiometer (HRIR) of the Nimbus II satellite 
while the third spectral region corresponds to the spectral region of one 
of the channels of the Medium Resolution Infrared Radiometer (MRIR) of Nimbus 
II and III, and to the window channels of the temperature humidity infrared 
radiometer (THIR) of Nimbus IV and V. 

The 6-8 ym region is located in a water vapor absorption band. It is 
in this spectral region that a number of sensors have been designated to 
obtain measurements from which atmospheric water vapor quantities are 
inferred. These radiometers have included the water vapor channel of the 
MRIR of Nimbus II and III and the water vapor channel of the THIR of Nimbus 
IV and V. 

The last spectral region covers a carbon dioxide band. This region has 
been used to obtain measurements of lower stratospheric water vapor as in the 
MRIR of Nimbus II. More recently, it has been used to obtain vertical profiles 
of atmospheric temperature. No analysis of this spectral band was performed here. 

3.3 Computations 

3.3.1 Procedure 

The computation program developed by Dr. David Anding was modified 
and used to evaluate the radiation transfer process through a cloudless non- 
scattering atmosphere. Details of the program may be found in the referenced 
document. Modifications to the program are given in Appendix B. Basically, 
the program uses a library of absorption coefficients for the various gaseous 
species and, with specifications of atmospheric profiles of temperature, absorber 
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concentration, and pressure-height, it uses a layered geometry model of the 
atmosphere to compute a number of radiation quantities at fixed wave- 
lengths for specified wavelength intervals. Since the absorption co- 
efficients are based on empirical data, the spacing of the wavelengths 
at which computations are performed changes with the wavelength region. 
Tables 3-1 and 3-2 show two examples of this spacing. The first is a 
print-out of the spectral radiance and atmospheric transmissivity com- 
puted for the 4 um water vapor window region using a specific atmo- 
sphere. The second table. Table 3-2, shows similar results for the 
10-12 ym window region. It is clear that the computations are made at 
much closer spectral spacing in the 4 ym region than in the 10-12 ym 
region. 


3.3.2 Parameters 


The Anding program was used to compute directly the following terms: 

_2 

Spectral radiance, watts/cm sr-ym 

Blackbody temperature, T (°K) 


Atmospheric transmissivity, x^(A) 


and, in the case of the absorption band regions, computations were also 
made, at selected wavelengths, of 


3x (z) 

Weighting function: * 

o Z 


(km -1 ) 


This term describes the relative contribution to the received radiance 
made by each level of the atmosphere. Figure 3-1 shows two examples of 
the weighting function computed for a wavelength in the 6.3 ym water vapor 
absorption band. In Figure 3- la, the computations were made for a very dry 
atmosphere. The weighting function shows that the atmosphere is nearly 
transparent with major contributions being made by and near the surface. 

The weighting function shown in Figure 3- lb was computed for a relatively 
moist atmosphere. The peak of the weighting, i.e., the level of maximum 
contribution, is shown located in the middle troposphere. It is clear 
therefore that the weighting function is useful in an analysis of the source of 
the radiant energy received by a satellite radiometer. 
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Table 3-1 

Example of Computer Printout for the 4 pm Region 


ATMOSPHERE ■ REGION b JAN SFNSQR ALTITUDE * 80000, SURFACE TEMPERATURE ■ 2T1.3* 

RUN NUMBER 10 HAS 221 WAVELENGTHS 


LAMDA 

DELTA 

RADIANCE 

TRANSMISSION 

3,83i 

0,0075 

i,i5«<Ma8P-05 

0.93529 

3 , 039 

0,0075 

l,3«67i37F-05 

0.9500? 

3,646 

0,0075 

l,a0B5B07r-05 

0,v««95 

3.654 

0,0075 

l,«290687F-05 

0.93796 

3,661 

0,0070 

l,«53876«F-05 

0,94224 

3,666 

0,0075 

t,«792082F-0S 

0.9^696 

3,676 

0,0075 

|,50<tnoi<lE-05 

0.94336 

3,663 

0,0075 

t.5?«?«?6F=C5 

O. 93717 

3.69J 

0,0080 

1.6525700F-05 

0,93994 

3,69V 

0,0075 

1.579742BF-05 

0.94O31 

3,906 

0,0075 

1,6O1R50BE-O5 

0,93667 

3,914 

o,oo«o 

1.6277907F-05 

0,93378 

3,922 

0,0075 

1.65310 W-05 

0.92904 

3,929 

0,0075 

1,6771 932F-05 

0,92940 

3,937 

0,0080 

l,70l94i25F-nb 

0.92411 

3,945 

0,0080 

1 ,7?7 a 79 7F-05 

0.92023 

3,953 

0,0080 

1, 752«55lE-05 

0,91482 

3,960 

0,0075 

l,77tt30«BE-05 

0,91084 

3,966 

0,0075 

l,7982195t-05 

0,90406 

3,976 

0,0080 

1.8222927F-05 

0.6982! 

3,984 

0,0080 

l,80««816F-05 

0.80955 

3,992 

0,0080 

l,e665523F-05 

0.88100 

4,000 

0,0080 

1.8862091F.-05 

0,86961 

4,008 

0 , 00«0 

1 , 9088 30 0F-05 

0.86?b9 

4,016 

0,0080 

1 .9299987F-05 

0.85375 

4,024 

0,0080 

1. 907761 1F-05 

0,64151 

4,032 

0,0080 

1.9607592E-05 

0.8230? 

4,040 

0,0080 

1.9821251F-05 

0,0156? 

4,049 

0,0085 

2.0098421F-05 

0.80573 

4,057 

0,0085 

2.0159336E-05 

0,7087? 

4,065 

0,0080 

2,02680?4F-O5 

0.77223 

4,073 

0,0080 

2,03337«0F-O5 

0,74777 

4,062 

0,0085 

2.0569685F-0S 

0.74122 

4,089 

0,0075 

2,072571 IE-05 

0,73305 

4,096 

0,0080 

2.0669178F-05 

0.71715 

4,107 

0 , 00*>0 

2.0987121F-OS 

0,69564 

4,115 

0,0085 

2.1015803F-05 

0.66933 

4,124 

0,0085 

2, 1502536F -05 

0.66757 

4,132 

0,0085 

2, 1737869F-05 

0.66146 

4,141 

0,0085 

2.2005406F-05 

0 , 67?9 1 

4,149 

0,0085 

2.2293636E-05 

0,67051 

4,156 

0,0085 

2.2592591F-05 

0,66665 

4,167 

0,0090 

2.2351058F-05 

0.63766 

4,175 

0,0065 

9,2«502O0F-06 

0.10063 

4,164 

0,0065 

2.09952a0E-06 

0.00001 

4,193 

0,0090 

1.6622139F-06 

0.0 

4,202 

0,0090 

1.6953536E-06 

0.0 

4,211 

0,0090 

1,7«00271E-06 

0.0 

4,219 

0,0085 

1.7757966F-06 

0.0 

4,226 

0,0085 

1.6156<l25F-06 

0,0 

4,237 

0,0090 

1 , 856 1 605F-06 

0,0 

4,246 

0,0090 

1.8973587F-06 

0,0 

4,255 

0,0090 

l,93919<l6E-06 

0.0 

4,264 

0,0090 

l,96l8«35F-n6 

0.0 

4,274 

0,0095 

2.0299876F-06 

0,0 

4,263 

0,0095 

2, 0740772F «06 

0,0 

4,292 

0,0090 

2.1108500F-06 

0,0 

4,301 

0,0090 

2, 1693773E-06 

0,0 
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Table 3-2 


Example of Computer Printout for the 10 uni Region 


ATMOSPHERE ■ REGION 5 JAN SENSOR ALTITUDE « 60000,' SURFACE TEMPERATURE ■ 2M,‘jO 


run number is has 122 wavelengths 


LAMOA 

DELTA 

RADIANCE 

TRANSMISSION 

9,070 

0,0700 

4,47387tlF-04 

0.97115 

9,130 

0.0750 

9.5065070E-O9 

0.97173 

9.220 

0,0900 

9.5522093E-09 

0.97205 

9.310 

0,0570 

9.5939959F-0U 

0.96537 

9,335 

0.0950 

9.6055629F-09 

0,96544 

9,398 

0,0630 

4, 0*S?954F -04 

0.*01 45 

9,463 

0,0990 

2.0*9?70lE-09 

0.16701 

9,494 

0,0310 

1,909261 IE-09 

0.13574 

9.526 

0,09*0 

2, 1002995E-09 

0,18164 

9,590 

0,0630 

2,8232979E-09 

0,36*69 

9,652 

0,0620 

2,166**937-09 

0.19585 

9,713 

0,0600 

2.39221 637. -0 9 

0.25062 

9,773 

0.0610 

2,979*373F-09 

0.26763 

9,839 

0,0600 

2, 6567653 F-OU 

0.31000 

9,893 

0,0590 

2.9069162E-09 

0.38332 

9,953 

0,0890 

3, 298*89 3E-09 

0.4928? 

10,070 

0,tl*0 

3,93339787-09 

0,66054 

10,190 

0,1100 

9,9072l*7F-09 

0 . *2580 

10.290 

0,1250 

9,90061 29E-09 

0,95*72 

10,990 

0,1300 

9, 924599 OF-OO 

0.96235 

10,550 

0,2300 

9,93699*OE-09 

0.95385 

10,900 

0,2500 

9,9511 157F-09 

0,99*64 

11,050 

0,1370 

9 , 99 33950F-09 

0.94117 

11,173 

0,1090 

9,92697297-09 

0.93314 

11.268 

0,0950 

9.9?896*3E-09 

0.93042 

11,369 

0,1160 

9,921*9567-09 

0.92606 

11,500 

0,1530 

9,916*3607.-09 

0.92073 

11,670 

0,1950 

9,91019917-09 

0.91623 

11,790 

0,1150 

9,69979527-09 

0.91461 

11.900 

0,1300 

4,679?82eE-09 

0.90*29 

12,050 

0,1500 

4.840*1 92F-04 

0 , *938* 

12.200 

0,1500 

4.7P9P805F-04 

0.87182 

12,350 

0,1500 

4 , 732*27 7F -04 

0.84749 

12.500 

0,0900 

4.6669063F-04 

0.81786 

12,530 

0,0300 

4,6601 P66F -04 

0.81064 

12.560 

0,0300 

4,64751 85 F -04 

0,80796 

12.590 

0,0 

4.5453501F-04 

0,71172 

12.630 

0,0 

4,51706*2E-04 

0,69446 

12.660 

0,0300 

4.48266747-04 

0,66198 

12.690 

0,0300 

4,6O77469E-04 

O. 80231 

12,720 

0,0350 

4, 5*81979 7. - 04 

0.79O56 

12.760 

0,0350 

4 , S70 086uE-04 

0.78257 

12,790 

0,0300 

4.5541720E-04 

0.77519 

12.620 

0,0300 

4.5373221E-04 

0.76681 

12,850 

0,0350 

4,51 7/829F-04 

0.75027 

12.890 

0,0350 

4.5011449E-04 

0.74813 

12,920 

0,0300 

4 , 4* 1221 5F -04 

0.73571 

12.950 

0,0350 

4 ,458241 2F-04 

0,72049 

12,990 

0,0350 

4.4349954E-04 

0,70811 

13.020 

0,0350 

4 , 4096960F-04 

0.69162 

13,050 

0.0350 

4.3839798F-04 

0.67682 


45 



•4 SI 



Figure 3-1 Examples of Weighting Functions Computed for the 6.3 pm Water Vapor Absorption Band 




Spectrally weighted values of the radiance were then computed using 
the results of Anding's programs. The spectral weighting consists of the 
following procedure: 

N = I * x N x di (3-2) 

A 

where is a spectral response function which essentially describes the 
relative efficiency of a real sensor to detect and measure incident radia- 
tion at various wavelengths. In this formulation, <p ^ has values between 0 
and 1. The spectral response function of the specific satellite radiometer 
channels selected for analysis are shown in Figure 3-2a to 3-2d. 

An additional parameter evaluated was the effective blackbody tempera- 
ture T,, (°K) corresponding to the computed effective radiance. These 

values were obtained by using the calibration charts for each of the radiom- 
eters as given in the Data User's Guides . Figure 3-3 shows examples of 
calibration curves for the HRIR of Nimbus II and the water vapor channel 
of the THIR for Nimbus IV. 

3.3.3 Atmospheres 

The atmospheres used in the data analysis include the following: 

1. Standard Atmospheres at 15°N, 30°N, 45°N, 60°N and 75°N for the 
months of January and July (with the exception of the 15°N model 
which is described by a single set of annual values) . The data 
for these models were obtained from the Handbook of Geophysics 
and Space Environments and the U.S. Standard Atmosphere Supple- 
ments , 1966. 

2. 4-D atmospheres from each of the 45 regions for the months of 
January and July. 

The results of these computations are presented in Section 4 of this 
report . 
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Relative Response % Percent Transmission 
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Wavelength fim 


Wavelength p m 


(a) Relative Spectral Response of 
3-4 ym Channel on NIMBUS II 


(b) Relative Spectral Response of 
10-11 ym Channel on NIMBUS II 



6.0 6.5 7.0 7.5 8.0 


Wavelength /im 

(c) Relative Spectral Response of 
6.7 ym Channel on NIMBUS IV 



11.5 ym Channel on NIMBUS IV 


Figure 3-2 


Examples of Effective Spectral Response Functions from NIMBUS II 
and NIMBUS IV (from NIMBUS II and IV User's Guide) 
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4. INFRARED ANALYSIS 


The computations presented in this section are meant simply to demon- 
strate the spatial variability of radiometric parameters introduced by the 
use of the MSFC 4-D Model atmospheres. The radiometric model by Anding 
used in this study can provide good simulation of the spectral radiance at 
any of the specified wavelengths. However, given the finite spectral band 
widths of specific radiometer systems, as defined by their effective spectral 
response functions, there is no clear and unambiguous approach to the simu- 
lation of the radiometric quantities actually measured by a real system. The 
approach adopted in this study (That of using the effective spectral response 
function to simulate the effective radiance, and using this computed value 
to derive the effective blackbody temperature from a calibration chart) is 
subject to a number of interpolation errors. Consequently, the computation 
results in Section 4 should not be used to construct absolute correction 
factors for the specific sensors discussed. More appropriately, the data 
should only be used to interpret the relative spatial variability of the 
radiometric quantities anticipated in simulation results using the MSFC 4-D 
Models. 

4.1 Introduction 

The MSFC 4-D atmospheric model showed marked differences in temperature, 
water vapor and pressure across its 45 homogeneous moisture regions. In fact, 
the analysis of Section 2 indicated significant climatic differences within 
latitudinal belts that are represented by only one atmosphere in the U. S . 
Standard Atmosphere Supplements . This section will look at the radiometric 
parameters derived from these same atmospheres to see if the significant 
atmospheric differences seen in Section 2 cause significant radiometric dif- 
ferences, or whether the standard atmospheres can provide an adequate repre- 
sentation of global atmospheric attenuation. 

Section 3 presented the radiometric model and computation techniques 
used to determine the blackbody temperatures presented in this section 
Nine radiometers will be examined: the Nimbus II 3-4 ym (HRIR) , 6-7 ym, 
and 10-11 ym (HRIR) channels and the Nimbus IV 6.7 ym (THIR) and 11.5 ym 
channels. The 3-4 ym channel, and the 10-12 ym channels are located in 
atmospheric windows; their values should be relatively independent of 
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atmospheric structure and composition, and highly dependent on surface 
temperature. The other two channels are located in the water vapor absorp- 
tion band; the variations in their values should show a high correlation with 
the atmospheric variations seen in the moisture regions. 

The blackbody temperatures will be presented for all the regions and 
standard atmospheres to permit a summary of global variability. They will 
then be examined for differences within the same latitudinal belts analyzed 
in Section 2 to determine the sensitivity of infrared transmission to atmo- 
spheric differences. Last, the radiometric values will be given for two 
typical suborbital tracks to allow an estimate of the actual variations 
encountered within one orbit. 

4.2 Overview of the Regional Radiometric Data 

Figure 4-1 presents the surface temperatures for the 45 regions for 
winter and for summer. It should be noted that the data for the Southern 
Hemisphere regions 37 to 45 for July appears with the Northern Hemisphere 
January regions and that the corresponding data for January is grouped with 
July. In the radiometric model, the earth was assumed to have an emissivity 
of 1; these temperatures are thus the effective blackbody temperatures at 
the earth's surface. The precipitable water values for these regions are 
the same as those discussed in Section 2; reference to those values, and 
to the plots of atmospheric structure, will be made throughout this section. 

4.2.1 The Window Channels 

The computed effective blackbody temperatures for the three "window" 
channels are presented in Figures 4-2 and 4i-3. It is immediately clear 
that some attenuation has occurred in all channels, even though they are 
located in atmospheric "windows". None of the temperatures given should 
be taken as precise values of effective blackbody temperature due to 
approximations in the calibration and conversion from spectral radiance 
to blackbody temperature (this is especially true for the 3-4 ym channel) ; 
however, the errors are consistent for each channel permitting comparisons 
between the regions and some evaluation of the atmospheric effects. 

The blackbody temperatures measured by 3-4 ym channel, and the two 
10-11 ym channels, are highly dependent on the surface temperatures beneath 
the atmosphere. A glance at the figures show that the radiometric tempera- 
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Figure 4-1 Region Su: 









tures generally increase with increasing surface temperature and decrease 
with decreasing surface temperature. However, all channels show the atmo- 
spheric contamination in that a complete correlation between the surface 
and radiometric temperatures does not exist; many changes are damped or 
accentuated by intervening water vapor. The two channels least attenuated - 
the Nimbus II 3-4 pm and the Nimbus IV 11.5 pm - are very close to the sur- 
face temperature in Regions 1, 2, 4, 5 and 37 during the winter. All of 

_2 

these regions have less than 0.5 g cm of precipitable water for these 

seasons. On the other hand, temperature differences of 3°K for the 3-4 pm 

channel, and from 4 to 6°K for the 11.5 pm channel for Regions 16 to 22 in 

_2 

January show the attenuation due to 2 to 4 g cm of precipitable water. 

The Nimbus II 10-11 pm channel shows a much greater sensitivity to 
water vapor than the two channels previously discussed. Even in the driest 
winter regions the blackbody temperatures are 5°K colder than the surface 
values; those regions with high moisture content often show a temperature 
drop of 15° to 20°K. The range of temperatures is decreased from the range 
of surface temperatures by 6°K in the winter and 8°K in the summer; the 
most marked changes are in the hot, humid regions where the maximum tempera- 
tures recorded are 289° rather than 300° K in the winter and 290° rather 
than 306°K in the summer. In July Regions 17 through 22 all have blackbody 

temperatures within 2°K of the one another; yet the surface temperatures 

_2 

increase 4°K and the precipitable water vapor goes from 4 to 5 g cm . On 

the other hand, in January the blackbody temperatures increase 7°K for a 

surface temperature increase of 8°K; in this season the precipitable water 

.2 

goes from 2 to 4 g cm . 

4.2.2 The Water Vapor Channels 

The two 6-7 pm channels presented in Figures 4-4 and 4-5 exhibit 
totally different characteristics from those found for the three previously 
discussed channels. At these wavelengths, the absorbtivity due to water 
vapor is so high that very large decreases in temperature occur, and the 
differences in surface temperature are almost totally obscured. The 
NIMBUS II channel runs 1-2°K colder than the NIMBUS IV with values between 
225°K and 255°K in the winter, and between 235°K and 255°K in the summer. 
Comparison of these blackbody temperatures with the precipitable water for 
the same regions reveals, as expected, a very high inverse correlation. To 
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Figure 4-4 "Water Vapor" Channel Blackbody Temperatures for th 
Winter Region 



Figure 4-5 "Water Vapor" Channel Blackbody Temperatures for the 
Summer Region 


obtain the fairly smooth blackbody temperature plots, some influence of the 
surface temperature is clearly evident; since the temperature differences 
increase from 15°K in the Arctic to 50°K in the tropics (and to 66°K in 
the monsoon season of Region 36), a compensatory increase of 39° to 50°K in 
surface temperature is necessary to maintain most of the blackbody tempera- 
tures in the 240° to 250° K temperature range. It is intuitively clear that 
the higher moisture values are usually associated with the higher surface 
temperatures; where this is not so (e.g.. Regions 33-35 July), significant 
changes in these blackbody temperatures do occur. 

4.3 Comparison with the Standard Atmospheres 

The blackbody temperatures computed for the standard atmospheres 
are given in Figures 4-6 and 4-7. Since only nine values represent all 
latitudinal belts for both seasons, both January and July appear together 
in the plot of each channel. For clarity, a plot of the standard atmosphere 
surface temperatures is given in Figure 4-7. 

The range of blackbody temperatures computed from the standard atmo- 
spheres duplicates well the range computed from the 4-D atmospheres; only 
the extremely hot values in the "window" channels corresponding to the 
slimmer deserts, and the extremely cold values in the "water vapor" channels 
corresponding to the Antarctic winter, are not reproduced. To permit 
comparisons between the two data sets, all blackbody temperatures computed 
from the region data were regrouped by latitudinal belts and appear in 
Figures 4-8 through 4-17. Since surface temperatures vary greatly within 
each zone, and since the differences between the surface temperatures and 
the blackbody temperatures are the values most representative of atmospheric 
attenuation. Figures 4-18 through 4-27 present the "window channel" tempera- 
ture differences applicable to each latitudinal belt. On the temperature 
difference plots, the value computed from the appropriate standard atmosphere 
appears as a constant across the graph; both sets of figures may be easily 
compared with the precipitable water plots given in Figures 2-14 through 2-23. 

4.3.1 The Arctic and Antarctic Zones 

The actual blackbody temperatures computed from the standard atmospheres 
for 75°N (Figures 4-8 and 4-9) correspond well to the values computed from 
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Figure 4-8 Blackbody Temperatures 60.75 
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Figure 4-10 Blackbody Temperatures 50.25° -60.75°N 
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Figure 4-11 Blackbody Temperatures 50.25 -60.75 S 








6104 



(Xo) aa i 



Figure 4-13 Blackbody Temperatures 37.5 -50.25 S 





Region 

Figure 4-14 Blackbody Temperatures 22.5° -37.5°N 


67 








Region 


Figure 4-16 Blackbody Temperatures 0° -22.5°N 
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Figure 4-18 
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Figure 4-19 "Window" Channel Blackbody Temperature Decrease 
60.75° -90°S 
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Figure 4-20 "Window" Channel Blackbody Temperature Decrease 
50.25° -60.75°N 
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Figure 4-23 "Window" Channel Blackbody Temperature Decrease 
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Figure 4-24 "Window" Channel Blackbody Temperature Decrease 
22.5° -37 . 5°N 
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Figure 4-26 "Window" Channel Blackbody Temperature Decrease 
0° -22.5°N 
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the Region 1 atmosphere and, accordingly, are too low for the values computed 
from the other regions north of 60.75°N. As expected from the atmospheric 
analysis in Section 2, the standard atmosphere gives radiometric values too 
cold for the Antarctic winter and too warm for the Antarctic summer (Regions 
37 and 38) . 

Figures 4-18 and 4-19 show that all regions north of 60.75°, with the 
exception of the warm, moist Region 3, show a temperature decrease of 0-2°K 
in the 3-4 ym channel. The 75°N standard atmosphere's value of 1°K thus 
represents this channel very well in these latitudes. Regions 37 and 38 
show a temperature decrease of 3°K; the standard atmosphere is slightly 
more in error here. 

The Nimbus II 10-11 ym channel shows temperature differences of 3 to 7°K; 
the value of 5° K computed for 75°N is a good median and accurately represents 
the Antarctic. The Nimbus IV 10-12 ym channel is not so well behaved; there 
the differences between the surface temperatures and the radiometric temper- 
atures range from 0 to 6°K, and the standard atmosphere's value of 4°K 
overestimates the attenuation for all but the ocean regions. Regions 6, 9 and 
38. Less attenuation is predicted for Region 4 than for the drier and colder 
Regions 1 and 2; this reflects, in part, the error in measuring small differ- 
ences in large numbers and, in part, the effects of atmospheric density 
and vertical temperature distribution. The standard atmosphere for July 
provides a good estimate of the attenuation expected for Regions 1, 2, 37, 
and 38 during the summer; however, the other regions in this zone are clearly 
subarctic in this season and are not well represented by the values for 75°N. 

4.3.2 The Subarctic Zone 

In the computed blackbody temperatures for the "window" channels 
(Figures 4-10, 4-11), the 60°N standard atmosphere resembles Regions 3 and 
6 during July, and no regions during January. As was seen in Figures 2-6, 
the standard atmosphere surface temperature is colder than any region sur- 
face temperature between 50.25 and 60.75° during the winter season. 

The predicted attenuation is far more representative of this latitudinal 
belt than are the predicted blackbody temperatures . No one region is 
accurately represented for all channels but many region values do not depart 
significantly from those predicted by the standard atmosphere (Figures 4-20, 
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4-21). This is especially true for the Nimbus IV 10-12 ym channel in July; 
the temperature loss indicated by the regions ranges from 4° to 7°K; the 
value of 6°K computed from the standard atmosphere is excellent. The same 
expected temperature decrease is given for the 3-4 pm channel for the 
same month; since 6 of the 9 regions have expected losses of 5 to 6°K, 
this value is also quite good. The standard atmosphere's 2°K temperature 
drop computed for January represents a good estimate as it is both a mean and 
a mode of the region values. For both seasons, the standard atmosphere 
slightly underestimates the attenuation found in ■ the NIMBUS II 10-llum 
channel; in January just the opposite occurs in the Nimbus IV 10-12 pm 
channel where attenuation is overestimated for all but Regions 3, 6, 27, 
and 28 . 

In keeping with the water vapor amounts shown in Figures 2-16 and 2-17, 
the standard atmosphere does not predict enough attenuation in the water 
vapor channels for any of the January regions but Regions 2 and 4; in July 
its values are much more typical. Since this is a transition zone, which 
shows an increase of 8°K in the January "water vapor" blackbody temperature 
from Regions 2 to 9, it is important that the region data be retained to 
provide adequate radiometric information. 

4.3.3 The Mid-Latitude Zone 

Most of the regions between 37.5 and 50.25°N have already been presented 
in the previous section. However, the additional regions match well with 
the regions already in evidence, and a glance at Figures 4-22 and 4-23 make it 
immediately clear that the attenuation predicted by the standard atmosphere 
for 45°N is representative for all channels but the NIMBUS IV 11.5 urn 
channel in July. There the estimated amount of attenuation is too large. 

The actual radiometric values are misleading because this, like the sub- 
arctic zone, is an area of highly variable surface temperatures (Figures 4-12 
and 4-13); a detailed breakdown of surface temperatures is necessary to obtain 
the blackbody temperatures. However, if an estimate of temperature loss 
is desired, the standard atmosphere would rarely be off by 3° K. 
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4.3.4 The Subtropics 


Because of the large climatic variability in the latitudes between 
22.5° and 37.5°, two-thirds of the 4-D moisture regions are represented 
in these zones. This leads, of course, to very large variations in the 
blackbody temperatures; the NIMBUS II 3-4 ym channel varies 45°K in 
temperature during the Northern Hemisphere winter (Figure 4-24). The 
summer values in the "window" channels are much more stable (the surface 
temperatures only vary by 25°K) but the high variability in moisture 
content, coupled with the varying surface temperature seen in both seasons 
(Figure 2-20), causes a 30°K range in the values of the water vapor channel. 
Normally, increased temperatures are matched with increased water vapor 
and thus to increased attenuation; in this latitudinal belt, the very hot, 
dry deserts are a clear exception to the rule. The Southern Hemisphere 
is much less variable at these latitudes; it is not, however, homogeneous 
(Figure 4-25). 

In light of the significant differences in temperature and moisture 
seen in these latitudinal belts, the plots of temperature decrease for the 
"window" channels (Figures 4-24 and 4-25) seem quite smooth. As expected, 
the winter values show a greater variability than the summer ones; yet 
the majority of the values lie within 3°K of the 30°N standard atmospheres. 

A notable exception is the Southern Hemisphere July NIMBUS IV 10-12 ym 
channel; here the standard atmosphere underestimates the attenuation in 
all regions. It is also slightly low for this channel during the northern 
hemisphere winter. Since the standard atmosphere overestimates the moisture 
for most of these regions, this result clearly illustrates the complex 
relationship between the atmosphere's structure and composition, and its 
attenuating effect on infrared radiation. The standard atmosphere also 
underestimates the attenuation in Regions 35 and 36 for all channels; 
these regions have a higher moisture content than that given by the 
standard atmosphere. 

4.3.5 The Tropics 

The tropics show less variability than the subtropics with a 30°K 
range of blackbody temperatures in the "window" channels indicating a greater 
similarity in surface values (Figures 4-16 and 4-17). The "water vapor" 
channels are more variable; this reflects the fact that the water vapor at 


83 


these latitudes varies more substantially than do the surface temperatures. 
Figures 4-26 and 4-27 show, however, that the annual standard atmosphere for 
15°N is not truly a representative value. It overestimates the attenuation 
for two of the three winter channels and underestimates it for two of the 
three summer channels. On the other hand, the variability in attenuation 
is so small that only the winter values for the Nimbus IV 10-12 ym channel, 
and the summer desert values for the Nimbus II 10-11 ym channel, consistently 
diverge from the standard atmosphere by more than 3°K. It would seem, 
though, that this zone might be better represented by separate atmospheres 
for summer and winter. 

4.4 Simulated Orbital Tracks 

4.4.1 Selection of Simulated Values 

The value of the atmospheric models in predicting the infrared tempera- 
tures sensed by satellites can be summarized by a simulation of typical 
orbital tracks. Since most satellites have a pole-to-pole orbit, two 
longitudes were selected for analysis, 90°W and 180°. The 90°W longitude 
approximates the center of the North American continent; in contrast the 
180° longitude intercepts land only at the tip of Siberia. This facilitates 
a comparison between oceanic and continental conditions while permitting 
a survey of differences (or similarities) across the latitudinal belts in 
both hemispheres. 

For this analysis, the regions appropriate to each longitude were 
selected, and their radiometric values plotted by latitude for each of the 
five channels (Nimbus II 3-4 ym, 6-7 ym, 10-11 ym. Nimbus IV 6-7 ym, 

10-12 ym) . These values appear as line drawings in Figures 4-28 through 
4-32. To show the blackbody temperatures that would be predicted by the 
standard atmospheres for the same tracks, the computed radiometric values 
were plotted as dots at the appropriate latitudes on the same figures. 

4.4.2 Analysis of Simulations 

The northern hemisphere is predominantly covered by land at 90°W 
longitude and is thus represented by the regions which most closely resembled 
the standard atmospheres (Section 2). This similarity reappears here; the 
standard atmosphere blackbody temperatures agree well with the 4-D tempera- 
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Figure 4-28 NIMBUS IV 10-12 ym Blackbody Temperatures at 90°W 
and 180° Longitude 
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Figure 4-30 NIMBUS II 10-11 ym Blackbody Temperatures at 
90°W and 180° Longitude 
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Figure 4-31 NIMBUS IV 6-7 ym Blackbody Temperatures at 90°W 
and 180° Longitude 
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Figure 4-32 NIMBUS II 6-7 pm Blackbody Temperatures at 90°W 
and 180° Longitude 
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tures for all latitudes except the tropics in the July Nimbus II pm channel. 

The radiometric values at this meridian also show the expected atmospheric 
differences; that is, the latitudes from 50° to 30°N show a temperature 
increase of 30°K in the "window" channel as the sensor travels south during 
January, and less than a 15°K increase during July. The blackbody temperatures 
for the tropics differ but little for the two seasons in the "window" channels; 
more of a change is visible in the "water vapor" channels emphasizing that the 
seasonal variability of water vapor here is greater than the seasonal vari- 
ability of temperature. 

The blackbody temperatures drawn from the two sets of atmospheres for 
180° agree within 2°K throughout the Northern Hemisphere. However, the 
standard atmospheres do seem slightly less representative at this longitude 
since the predicted temperatures are generally too warm for the July "window" 
channels, and too cold for the more northerly latitudes in the January "water 
vapor" channels. The only really significant difference between the two 
longitudes occurs in the latitudinal zone of 50° to 40°N; there the winter 
ocean blackbody temperatures run warmer than the land values, and the summer 
ocean temperatures run colder. This agrees well with the analysis of the 
temperature and moisture parameters done earlier. 

The standard atmospheres do not represent the southern hemisphere at 
either longitude except at 30 # S in July or, possibly, in the water vapor 
channels. During the winter a smaller range of temperatures is found than 
was seen in the northern hemisphere, and the transition from very cold to 
moderate warm temperatures is much more pronounced between 50° and 90°S 
than between 30° and 50°S. The summer range of temperatures is slightly 
greater than that indicated by the standard atmospheres but all temperatures 
are colder. This is partly due to the dominance of ocean regions on this 
half of the orbital subtrack; it is more due to the different climatic patterns 
of this hemisphere. 
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5. DISCUSSION OF RESULTS 


5.1 Comparison with Other Model Computations 

To demonstrate the validity of the computations made with Anding's model, 
the results for Supplementary Atmospheres were compared with computations made 
by other authors. Salomonson (1970) used the same Supplementary Atmospheres 
to calculate the expected effective blackbody temperature for the HRIR (with 
spectral response in the 3.4 pm region). The radiometric model used, however, 
was one developed by Kunde (1967). The results of his calculations are shown 
in Table 5-1, together with the corresponding value obtained in this study. 

The table shows that the present computations, with the exception of the 
tropical results, appear to be systematically colder than the results obtained 
by Salomonson. This might be due to a systematic error in the interpolation 
of the effective spectral response functions at specific wavelengths. However, 
the differences of 1°K and 2°K for the tropics are certainly within accept- 
able limits. Apart from agreement in absolute values, the more important 
aspect of the comparison in Table 5-1 is the trend in the data. In both 
sets of computations, the seasonal and regional variability is maintained. 

Another set of computations useful as comparison with the results pre- 
sented in this study are the ones recently reported by Stevenson and Miller 
(1972). In this study, preliminary computations using Kunde* s technique 
were performed of the expected values of AT corresponding to four meteoro- 
logical conditions. The results for the downward looking case are shown 
in Table 5-2. 

The importance of these results is that the computed temperature dif- 
ferences are within the limits shown in Section 4 of this study, and that 
the expected atmospheric effect is greater for the MRIR than for the HRIR, 
a result clearly brought out in the data shown in Section 4. Perhaps more 
important is the fact that the data shown illustrates that even diurnal 
effects can produce very significant differences in the expected radiometric 
values, a factor which must be considered in the simulation of such values. 
With the variance statistics built-in the MSFC 4-D data, these small scale 
effects can be effectively simulated by a properly designed Monte Carlo scheme. 
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TABLE 5-1 


COMPARISON OF SIMULATED EFFECTIVE BLACKBODY TEMPERATURES FOR 
THE HRIR USING SUPPLEMENTARY ATMOSPHERES 


Atmosphere 


1) Tropical 

2) Subtropical 
Summer 

3) Subtropical 
Winter 

4) Mid-Latitude 
Slimmer 

5) Mid-Latitude 
Winter 

6) Sub-Arctic 
Summer 

7) Sub-Arctic 
Winter (Cold) 

8) Arctic 
Summer 

9) Arctic 
Winter (Mean) 


Temperature 

Surface 

■ C°K) 


Effective Blackbody 
Temperature, T^ (°K) 


"Salomonson" t "Anding" 


299.65 

294.37 

296 

301.15 

295 

297 

287.15 

284 

283 

294.15 

290 

289 

272.59 

271 

270 

287.15 

283 

282 

257.28 

257 

: 

256 

278.15 

1 

275 

276 

249.22 

249 

248 


TABLE 5-2 

AT« (T - T bb ) VALUES DERIVED FROM STEVENSON 8 MILLER (1972) 


Time 

,<2pr-cm of . 
"Dry"( water va por^ 

r .>2pr-cm of 
"Wet" Skater vapor' 

Sensor 

0000 Local 

1200 Local 

0000 Local 

1200 Local 

HRIR 

2. 5°K 

2.5°K 

3.3°K 

3.2°K 

10-12 ym MRIR 

4.5°K 

5.0°K 

6.7°K 

5.8°K 
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5 . 2 Summary 


From the detailed analysis of results presented in Sections 2 and 4, the 
following points are noted: 

1. Along latitudinal belts, the effects of "continental ity", i.e., 
the presence of land and water masses, can introduce large varia- 
tions in moisture content within a latitudinal belt. 

2. For each of the latitudinal belts studied, the large spatial dif- 
ferences between regions are found during the summer months, where 

_2 

the range of precipitable water can be in excess of 2 g cm . For 

_2 

the winter the range is generally less than 1 g cm 

3. As expected, the most pronounced inhomogeneities in precipitable 

water are found in the July data (both hemispheres) for the 

latitude belt between the equator and 22.5°N(S). Here the pre- 

cipitatable water ranges from a low value of about 1.5g cm -2 to 

_2 

a high value of about 6.6 g cm . 

4. For the latitude regions poleward of 37.5°, the values of precipi- 
table water computed from the standard atmosphere appear to be 
representative of the regional mean for the corresponding latitude 
belts. Equatorward from 37.5°, the supplementary atmosphere 
appears to be systematically too high to be representative of the 
mean for the corresponding latitude belt. 

5. These differences in moisture between the supplementary atmo- 
spheres and the 4-D atmosphere are carried through to the radio- 
metric quantities, but to a lesser degree. While the range of 
precipitable water for a given latitude belt often exceeds 50% 
of the mean value, the corresponding radiometric temperature 
seldom exceed 1/30 of the mean values of the effective blackbody 
temperatures for the latitude belt. 

5.3 Conclusions 

The analysis in this study shows that the supplementary atmospheres as 
given in the Handbook of Geophysics are probably more representative of the 
temperature structure of the regimes for which they are designed than they 
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are of the moisture structure. The most immediate difficulty using the 
Supplementary Atmospheres is the high differentiation in the moisture field 
between land and water, a characteristic not included in the supplementary 
atmospheres. This differentiation is most pronounced in the tropical belts 
with their wet equatorial forests and extremely dry deserts. 

These regional differences in moisture field are reflected, although to 
a lesser degree, in the differences in simulated effective blackbody tempera- 
tures. The use of the supplementary atmospheres alone, especially in the mid- 
ana low-latitude regions, fails to bring out the significant patterns in 
expected effective blackbody temperatures in a given latitude belt. It is 
clear, therefore, that the limitations of atmospheric effect correction 
tables such as those presented by Salomonson (1970) and others must be 
recognized. For more precise corrections (to be used for local regions), 
detailed correction tables should be computed using the 4-D Model corre- 
sponding to that region. Thus, the correction factors can be made to reflect 
the true variability of the real atmosphere. 
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APPENDIX A 


A DISCUSSION OF WATER VAPOR PARAMETERS 
AND THEIR RELATIONSHIPS 

A. 1 Mathematical Relationships 

Considerable confusion has arisen due to the varying definitions and 

expressions of water vapor. The form used in the MFSC 4-D model is absolute 

humidity, p , also known as water vapor density, and is expressed in units 
-3 w 

of g m . A more commonly used unit is mixing ratio, w, a dimensionless 
term usually expressed in g kg . To facilitate conversion, the relation- 
ship between the two terms is given here as : 

W = (p u /p d ) * 10 3 (A-l) 

where 

W is the mixing ratio, g kg -1 , 

is the water vapor density, g m" , 

_3 

P d is the density of the dry air, g m 

In the 4-D model, p d is referred to simply as the atmospheric density; 
thus the mean mixing ratio is simply the ratio of the water vapor density 
to the atmospheric density. 

The dew point can be found from the mixing ratio graphically by plotting 
on the standard pseudo-adiabatic charts. When the atmospheric density is 

_3 

less than 1000 g m , and the absolute humidity is small, the dew point can 
also be approximated by using the absolute humidity value in place of the 
mixing ratio. This does, however, cause grave errors near the surface. 

_3 

There the density is greater than 1000 g m and the moisture is relatively 
high; use of the absolute humidity in such circumstances can overestimate 
the dew point by 3-4°K and lead to an erroneous indication of supersaturation 
or atmospheric instability. 

Computational means are necessary to obtain the dew point directly from 
the absolute humidity value. The equation requires knowledge of the atmo- 
spheric temperature, and is: 
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(A-2) 


T d * 237.3/(7.5 - l°g 10 ^1324.037^ ^ 

where is the dew point temperature, °C, 

T is the atmospheric temperature, °K, 

P oj is the absolute humidity, g m~ 3 

If it is desired to express the water vapor of the 4-D model in dew point 
temperatures, this is the most direct and accurate approach. 

A. 2 Example of Error in Dew Point Approximation 

Region 11, January, was selected to illustrate the dew point error 
caused by using the absolute humidity value for the mixing ratio value. 

The following mean values were given for the surface: 

Temperature 284.65 °K 

Pressure 1014.24 mb 

_3 

Absolute Humidity 8.38 g m 

Density 1241.33 g nf 3 

The last two terms gave a mixing ratio of 6.75 g kg * . When the 
absolute humidity and the mixing ratio were plotted on a pseudoadiabatic chart, 
they gave the following results: 

Dew Point (from absolute humidity) 284.8 °K 
Dew Point (from mixing ratio) 281.6 °K 

The use of the absolute humidity value not only overestimates the dew 
point by more than 3°K; it indicates that this region is supersaturated. 
Although supersaturation is not impossible for an ocean region, an average 
state of supersaturation for a five-year period is totally unreasonable. 
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APPENDIX B 


MODIFICATIONS TO THE PROGRAMS OF 
DR. DAVID ANDING 

Two programs were received from Dr. David Anding, CATM and RADST. 

The details of these programs are given in the report by Dr. Anding and 
are generally appropriate here, as no changes were made to the computations 
or theories given. It was necessary, however, to change the control cards 
to make the programs run on the 360/75 available to ERT. 

The major modifications were in the input and output portions of the 
programs. The program to set tip the atmospheric models needed for the 
radiation calculations, CATM, acquired new keywords and subroutines which 
input data from the 4-D region coefficients or from the ERT-NWRC disk files. 
These routines greatly increased the flexibility of the program and signi- 
ficantly minimized the time and effort needed to input these data sets. 

The program to compute infrared attenuation due to atmospheric absorption, 
RADST, was altered to print out vertical profiles of the weighting function 
for selected wavelengths, and to write on disk (or tape) values of total 
transmission as well as spectral radiance. 

All other changes were minor. For example, CATM was changed to inter- 
polate starting at the given surface height, if the initial level requested 
was below the surface height. This was important to retain the surface 
layer; previously, the output atmosphere would commence at the next 
highest level, i.e. 500 m, often losing a significant amount of moisture. 
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APPENDIX C 

AN EVALUATION ON THE EFFECTS OF CLOUDS ON STATISTICS 
OF PRECIPITABLE WATER 

C.l Introduction 

The statistics of atmospheric parameters which make up the 4-D Model 
Atmospheres were compiled independent of cloud cover conditions. For most 
climatological applications such cloud-cover-independent statistics are 
appropriate. Yet there exists a class of applications in which sub- 
sets of these statistics, those describing clear sky conditions only, 
are useful. One important area where such statistics are useful is the 
simulation of earth viewing missions by aircraft or spacecraft-borne visible 
and infrared sensors. In such an application the presence of clouds automa- 
tically precludes further consideration of atmospheric attenuation effects. 

The statistics of atmospheric temperature and moisture become important only 
when there are no clouds. 

Intuitively it can be argued that there should be differences in the 
structure of the atmosphere, i.e., its vertical temperature and moisture 
distributions, between clear sky and cloudy conditions. This is especially 
true in the case of the moisture field where the presence of clouds is an indi- 
cation of saturation at certain levels. With the limited variance in the 
temperature fields for a given climatic regime, the saturation in turn is 
an indication of an atmosphere more moist (in terms of absolute moisture 
parameters such as precipitable water or absolute humidity) than one with- 
out clouds. The means of precipitable water computed for all days would 
in fact overestimate the mean of precipitable water for clear sky days only. 

From the analysis presented in the report, it is clear that while 
the effective blackbody temperature or atmospheric transmissivity is not 
a simple function of any of the moisture parameters, the precipitable 
water determines to a large extent the degree of atmospheric interference 
in a cloudless atmosphere. Given this consideration, it becomes important 
to evaluate the degree to which the 4-D model moisture statistics are 
overestimates of those for clear sky days. In the following paragraphs, 
the results of some analyses directed towards this problem are presented. 

These results are compared with some findings by other investigators. 
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C.2 Analysis 


Data from three stations were selected for analysis. San Diego (03131) 

and San Francisco (23230) from moisture Region 18, and Tucson (23160) from 

moisture Region 30 were selected as being representative of regions with 

"moderate to high" and "low" annual mean moisture content. The period of 

record used in the analysis was the two year period between 1 January 1965 

and 31 December 1966. Daily values of precipitable water at each of the 

stations were computed using the rawinsonde data on magnetic tape. These 

values were then grouped by season (with winter being defined as December, 

January and February) and stratified by cloud cover conditions. To do 

this, use was made of the daily surface observation reports at 12Z made at 

each of the stations. Seasonal and annual statistics of the means and 

variances of precipitable water were then compiled for the following cloud 

cover conditions, some of which are overlapping (in tenths): 0-1, 2-3, 

4-9, 10, 0-3, 4-10, 2-10, 0-10, and "Precipitation". Results of these 

calculations are shown in Tables C-l, C-2, and C-3. Using these results, 

computations of the seasonal and annual values of the following ratios 

were computed: R. * W(0-l)/W(0-10) and R 9 = cr (0-1 )/o (0-10) where 

X fc w w 

W(0-1) and o (0-1) are the mean and variance of precipitable water for 
w 

days with cloud cover less than one tenth; and if (0-10) and o w (0-10) are 
the corresponding values computed for all days. These ratios are shown 
in Table C-4. 

C.3 Discussion 

Statistics of seasonal mean values of precipitable water for all 
sky conditions (cloud cover category 0-10) for the three stations show 
characteristically high values during the summer and low values in winter. 
The lowest seasonal values of precipitable water are found in Tucson 
during the winter months. The ratio of clear sky days to all days is 
approximately 1/2 for Tucson while the same ratio for San Diego and San 
Francisco is less than 1/3, reflecting the differences in climatic 
regimes . 

The annual and seasonal mean values of the ratios R^ and R 2 shown in 
Table C-4 do not indicate any significant differences between individual 
stations or climatic regimes. R^, the ratio corresponding to the seasonal 
mean values of precipitable water, ranges from a low value of 0.705 (Tucson/ 
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TABLE C-4 

SEASONAL AND ANNUAL VALUES OF Rj AND R 2 
(Rj = W(0-l)/W(0-10); R 2 * a w (0-l)/a w (0-10)) 



San Diego 

San Francisco 

Tucson 

Season 

R 1 

R 2 

R i 

R 2 

R i 

R 2 

Winter 

0.770 

0.405 

0.770 

0.657 

0.800 

0.448 

Spring 

0.892 

1.538 

0.825 

0.858 

0.929 

0.667 

Summer 

0.956 

0.748 

0.923 

0.457 

0.705 

0.608 

Fall 

0.893 

0.969 

0.796 

0.915 

0.831 

0.469 

Annual 

0.846 

0.948 

0.794 

0.846 

0.745 

0.403 
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Summer) to a high of 0.929 (Tucson/Spring) with most of the values being 
between 0.8 and 0.9. The values of I^, the ratio for the variance, show 
greater variability. It has a low value of 0.405 (San Diego/Winter) and 
an anomalously high value 1.53 (San Diego/Spring) . 

While it is not possible to draw definitive conclusions from results 
of the limited analysis presented here, it is nonetheless possible to tenta- 
tively conclude that: 

1. Clear sky days have a mean value precipitable water which is 
less than that obtained for all days. 

2. A "practical" factor to be used to adjust the 4-D mean precipi- 
table water data to obtain the subset of statistics for only 
clear sky days is most probably between 0.8 and 0.9. 

3. The variance of the mean precipitable water for cloudless days 
is (with minor exceptions) less than the variance for all days. 

4. The variability of the variance is sufficiently large as to make it 
difficult to arrive at an estimated representative factor 

to use in adjusting the 4-D data. It suffices to say that 
the 4-D values can greatly overestimate the variance for 
clear sky cases. 

It is perhaps significant to note that the findings in the present 
study regarding the means of precipitable water are supported by other 
investigations. Fritz (1949) suggested that the subset of annual 
mean of precipitable water at continental stations for clear sky days 
may be obtained from the set of annual means of the same parameter at 
the corresponding stations by multiplying each value by 0.85. Humphreys 
(1911), in an earlier study, arrived at a value for the factor of 0.89 
using European Stations. These values are certainly in consonance with 
the findings of this study. 
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